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Abstract 

Local scale invariance for lattice models is studied using new realizations of the Schro- 
dinger algebra. The two-point function is calculated and it turns out that the result can be 
reproduced from exact two-point correlation functions evaluated in the stationary state of 
several simple stochastic models. 



1 Introduction 

Among the most important concepts of physics is that of symmetries. Particularly in the context 
of critical phenomena, scale invariance and its consequences have been thoroughly studied for a 
long time. More recently, for isotropic, equilibrium critical points in two dimensions, conformal 
invariance has allowed for significant progress in the classification of the universality classes and 
the elucidation of their properties, see |1[ and references therein. 

Much less is known if the system under study is either strongly anisotropic or time-dependent, 
or both. Such systems will at criticality show an anisotropic scaling of, say, a correlation function 
C(r, t) depending on a "space" coordinate r (we take one space dimension throughout) and "time" 
coordinate t 



where 6 is the anisotropy or dynamical exponent and x is a scaling dimension. One may ask 
whether this strongly anisotropic scale invariance can be sensibly generalized to include space- 
time dependent rescaling functions A(r, t) instead of merely taking A to be a constant. Indeed, 
at least for systems for which the equations of motion reduce in the continuum limit to the free 
Schrodinger or diffusion equation (and thus have 6 = 2), this seems to be possible. The group of 
scaling transformations is then the Schrodinger group, which is the maximal kinematic invariance 
group of the free Schrodinger equation ||. The Schrodinger group can be obtained as the non- 
relativistic limit of the conformal group ||. 

The Lie algebra of the Schrodinger group in d = 1 space dimension is spanned from the set 
{X_i, X , Xi, Y_i, Yi, M } and the non-vanishing commutation relations are 



C(Xr, X 9 t) = \- 2x / e C(r } t) 




2 




(n — m)X, 



n+m 






2 J 







Note that M commutes with the entire algebra and the Casimir operator is M 



Q = (4M X - 2{Y_i , Yi }) 2 - 2{2M X_! -Y\, 2M X 1 - Y?} 

2 2 2 2 



(1.3) 
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We shall denote the eigenvalues of M by — M. and of Q by Q. The projective unitary irreducible 
representations (puir) of fli~2] ) are classified ||. There are four different types of PUlRs, but only 
one of them allows M to be realized non-trivially. Of these, the PUlRs which have Q ^ give 
rise to an infinite set of internal quantum numbers (spin does not arise in d = 1). We consider 
the PUIR with Q = which is realized on (scalar) wave functions ip{r,t) which are solutions of 
the Schrodinger equation {id t + l/(2Ai)d^)ip(r, t) = 0. It takes the form ip — > Uip, where [|J 

tt i / \ , 1-1/2 f i-M- f 1 2 St ~ P i i \ W 

Uip{r,t) = — yt + a\ 1 exp < — r H v — vr + (— 71 + a)av > 

{ -jt + a \2 2 / J 

, /r — (— ■yt + a)a — (5t — 3)av St — 3 \ . 
• ^ — ! 1 — — , — (1-4) 

with a5 — /?7 = 1 and a, /3, 7, 5, f and a are constants. This realization was found by Niederer 0] 
as giving the maximal kinematic invariance group of the free Schrodinger equation and, as shown 
by Hagen ||, also arises when discussing the Schrodinger invariance of non-relativistic free field 
theory. 

The Lie algebra generators are readily found from ( |1.4|) and read 

x n = - — rv<9 r - n{n + l) Mt n ~ l r 2 

2 4 

M = -M (1.5) 

X-i and Y_i generate translations in time and space, respectively, X generates global scale trans- 
formation, Yi is the generator of Galilei transformations and Xi generates the special Schrodinger 

2 

transformations. In order to implement Galilei invariance of the free Schrodinger equation, the 
scaling fields 0(r, t) of the theory will under a Galilei transformation pick up a complex phase 
characterized by the mass A4 and thus has necessarily to be complex ||, as is well known. Thus 
the field 4>(r,t) is characterized by its scaling dimension x and its mass A4, while its complex 
conjugate <p*(r,t) is characterized by the pair (x, —A4). This structure is to be kept when going, 
via the continuation A4 — > (2iD)^ 1 , to a free diffusion equation with diffusion constant D. The 
requirement of covariance of the two-point function under the given Schrodinger transformation 
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determines the two-point function (written here in Euclidean time) completely One finds 

< (j)i(r 1 ,t 1 )(f)* 2 (r 2 ,t 2 ) >= ®o$M 1 ,M2 S x 1 ,x 2 (ti -t 2 )~ Xl exp 

where <3>o is a normalization constant. Note that the exponential factor and the Bargmann su- 
perselection rule M.\ = M. 2 [0 follow from covariance under the Galilei sub-group spanned by 
{X_!,y ± i,Mo} alone. Several simple models can be shown to reproduce this result. Three-point 
functions and plane surfaces have been considered as well 0. 

These results for the correlation functions do depend not so much on the algebraic structure 
but rather on the particular realization of the Schrodinger algebra used. For example, one could 
give another realization by using the space of wave functions of the harmonic oscillator rather 
than free particles as a carrier space 0, [|. Here, we want to consider the consequences of another 
realization of ( |1.2|) . It differs from the free particle realization defined above in that the carrier 
space is the space of wave functions which are periodic in momentum space. We shall define this 
new realization in the next section and shall use it to obtain the two-point function by demanding 
covariance. Exact results for two-point functions directly obtained from some simple lattice models 
coincide with the expressions derived in section 2 and show that this realization can be viewed as 
evidence for a Schrodinger symmetry of lattice models. 

2 Lattice realization of the Schrodinger algebra 

We begin our discussion by noting that generators which satisfy the Schrodinger algebra (|1.2j) can 
be constructed from the generators {P, X, E, T, M} of a much simpler algebra A the generators 
of which satisfy the relations 

[P,X] = [E,T] = 1 (2.1) 

and 

[A,B)=0 (2.2) 

for all other pairs of generators A, B G {P, X, E, T, M}. The generators X 0j ±i, Y ± i, M defined as 
shown in table 1 satisfy indeed the commutation relations (|1.2|) . Therefore any representation of 
the algebra A gives rise to a representation of the Schrodinger algebra. Clearly, P = d r , X = r, 



Mi (n - rg) 
2 U-U 
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Table 1: Realization of the Schrodinger algebra in terms of generators of A. 



E = d t , T = t and M = M. (where r, t, M. are numbers) satisfy relations (|2.1|) and Q2.2|) and give 
rise to the "continuum" realization (|1.5|) of the Schrodinger algebra (see also table 2). 

The point we wish to make is that instead of choosing P = d r , one may represent P by the 
symmetric difference operator 

P = - sinh (^d r ) (2.3) 



a \2 

(which is defined by its expansion in a power series). In order to preserve relations ( [2.1| ) one may 
now choose 

The generators X ^±i,Y ± i, M$ defined in table 1, but now given by the expressions shown in 
table 2 still satisfy the commutation relations (|i~2"|). The various operators are understood to be 
defined in terms of a power series expansion. 

Some intuitive understanding of the role of the parameter a comes from considering the gen- 
erator Y_i of "translations" 

2 

Y_if(r) = ~ ( f(r + a/2) - f{r - a/2)) . (2.5) 

which, if we interpret a as a lattice constant, may be understood as a discretized symmetric lattice 
derivative operator. 
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-d t 

t 2 dt J* /9fl ,rsinh(a/2<9 r ) " ( u \ 9a ,r) 2 

1 acosh(a/2 a r ) \ ' ' ' 2 \ cosh(a/2 a r ) / 

— - sinh(a/2 ci.) 
2t sinh(a/2a r ) ./j, s ,r 

a V / ' / cosh(a/2 a r ) 
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Table 2: The "continuum" and "lattice" realizations of the Schrodinger algebra. 



The Casimir operator of the Galilei subalgebra has the abstract form C = ME — P 2 /2 and is 
therefore in the lattice realization given by 

C = Md t - 4 (cosh (a^) - 1) . (2.6) 
a 2 

Thus C can be viewed as the continuum or lattice Schrodinger operator in the "continuum" or 
"lattice" realizations of table 2, respectively. We emphasize that the generators given in table 2 
do not generate finite transformations but may still be viewed as generators of infinitesimal trans- 
formations of some Lie group. In this paper, however, we do not want to address the question of 
exponential maps of such generators. 

We now turn to the two-point function 

F = F(ri,r 2 ;ti,t 2 ) =< Mn, *i)^(r 2 , h) > (2.7) 

and study the consequences of covariance of F under the "lattice" realization of ( |1.2|) . We will 
use the short-hand 

In slight abuse of language we shall furthermore refer to the action of the operators Y m and X n as 
defined in table 2 as translation operator, Galilei operator etc. Two-particle operators are defined 
as the sum of two one-particle operators acting on the pair of coordinates (ri,ti) and (r 2 , i 2 ) 
respectively. 



Demanding time translation invar iance implies F = F(r 1 ,r 2 ; r) where r = t\ — t 2 . Invariance 
under M gives the Bargmann superselection rule M.\ — Ai 2 — -M.. Translation invariance, 
generated by Y_i, requires 



sinh ( i^D^j + sinh Q-A 



F(r 1 ,r 2 ;r) = . 



(2.9) 



This equation can be solved by introducing the coordinates R = r 1 + r 2 and p = r\ — r 2 and 
writing F = f(R, p; r). We get 



cosh (^d p J ■ sinh (^d R ) /( /?. r) = . 



This has the solutions 



f(R+%,r,T)-f(R-±r,T 



o 



(2.10) 



(2.11) 



and f(R,p + a/2;r) = —f(R,p — a/2;r). Since we would like to recover in the limit a — > 
the familiar continuum solution, we have to choose ( [2.11| ). The second choice would lead to 
F(ri, r 2 \ t\, t 2 ) = for a — > and we do not follow this possibility any further. 
Next, we consider Galilei invariance 



I -U sinh (£a) + tiMi — t~t~~ 
~l [a \2 ) cosh(a 



/2A 



-Ti \F = 



(2.12) 



where ei = — e 2 = 1. Using eq. ( [2.1 1| ) and expanding the cosh, we find 

1 



cosh(a/2<9 p 



sinh (ad p ) + M.p 



/(i2 + |,p;r) =0 . 



(2.13) 



Consider now the equation (cosh(a/2<9 p )) -1 /i(p) = 0. We write h in the form h(p) = Y^kLo hk{p)a k - 
Insertion then shows, after expanding the cosh, that order by order hk(p) = for all k. We are 
thus left with 



T 



sinh (ad p ) + Mp 



f[R+-,p;r 







(2.14) 



This can be solved via Fourier transformation with the result 



f[R+-,p;r 



g ( R+ -, T )— dkeir/Ma^^ikip/a) 

V 2 / 2n J--* 



--: g[R+^,T)$(p/a,T/Ma 2 ) 



(2.15) 



If we interpret the indeed to denote the sites of a lattice, p/a = n is an integer and a, rj .Ma 2 ) 
I n (r I Ada 2 ) , where J n is a modified Bessel function. We also note the following properties of $ 

1 



<9 r $(p/a, t/Mcl 2 
tanh (^<9 P J (p$(p/a, r/Ma 2 ) 



— — - cosh( ad p )$( p/a, r/Ma 2 ) 
M.a l 

-J^sinh 2 gd p )$(p/a,r/.Ma 2 ) 



(2.16) 



We now consider the scale transformations generated by X . Scale invariance requires that 

(2.17) 



2 f 1 1 

^ [ a cosh (a 



— r* sinh (- a") + -a* 1 F = 
/2A) V2 7 2 J 



where are the scaling dimensions of the fields 0j. Following the same lines as before and using 
that [p, f(ad p )] = —af'(ad p ) for any (differentiable) function /, one obtains, with x = {x\ + x 2 )/2 



* + ±-.(l*),-§ + , 



g{R,T)<&(p/a,T/Ma 2 ) = 



(2.18) 



Using (|2.16|) , this can be reduced to 



rd T + — T — + x - - 
Ma 2 2 



g(R,r) = 



with the solution 



g(R,r) = g (R)T 



-(z-1/2) 



exp 



(2.19) 



^2.20) 



Ma 2 , 

where go{R) satisfies go(R + a/2) = go{R — a/2). Finally, we demand invariance under the special 
Schrodinger transformation generated by X\\ 

2 r . . / , \ 2 



Ti sinh(a/2A) — 



i=l 



a cosh(a/2A) 



2 \cosh(a/2A 



- Xi ti } F = 



(2.21) 



A straightforward but tedious calculation shows that the solution found so far is consistent, pro- 
vided only that 

Xl = x 2 = x (2.22) 

The fields 0i and <fr 2 must thus have the same scaling dimension x the value of which is not 
determined by Schrodinger invariance. Summarising, the final result is 



F(n,t) =F(r 1 ,r 2 ;t 1 ,t 2 ) = g S Ml M2^ 1 ,x 2 t 1/2 Xl e *J n (t) 



(2.23) 



with 

t= t -^- and n= r -^± (2.24) 

Mia 2 a 

where we assumed that both r\ and r 2 are integer multiples of the "lattice" constant a. As in con- 
tinous space, Schrodinger invariance on the lattice determines the two-point function completely 
up to the non- universal amplitude go, the non- universal constant mass A4 and the universal crit- 
ical exponent x. Note that if we require 4>i an d 4>2 to be locally conserved, i.e. d t J2 n F{n, t) = 0, 
then one has x — 1/2. If r± or r 2 are not integer multiples of a, the Bessel function I n (t) has to be 
replaced by the scaling function ( |2.15 ). We can consider g Q as a normalization constant if both r\ 



and r 2 are integer multiples of a. 

3 Exact two-particle correlations in lattice models 

So far we have merely assumed the existence of physical systems in which lattice Schrodinger 
invariance holds and studied some of its consequences. Clearly it is desirable to have an example, 
even a very simple one, where correlation functions as predicted in the last section would be found. 
Such an example is many-particle lattice diffusion. 

3.1 Symmetric, non-exclusive lattice diffusion 

Consider a system of arbitrarily many particles hopping stochastically to their nearest neigh- 
bouring sites. We choose the probability of hopping from site i to site % ± 1 proportional to the 
occupation number of the site i. For this stochastic process one may either write a master equa- 
tion or, equivalently, a quantum Hamiltonian defining the time evolution of the system [|K| |ll], |12 



In this mapping the state of the system at some time t is given by (stochastic) occupation num- 
bers n = {rii} where i is the number of a lattice site. Their dynamics are given by a Schrodinger 
equation for the probability distribution F(n,t) with a quantum Hamiltonian H 

d t \F(n,t)) = -H\F(n,t)} (3.1) 

which for our very simple toy process is given by 

H = \ H( a I+i _ a !)0i+i - o-i) ■ ( 3 - 2 ) 



Qj j CLj 



8ij. Each state in the 



The operators a\ and satisfy bosonic commutation relations 
system is represented by a state vector \n) which together with the transposed vectors {n\ form 
an orthonormal basis. The creation operator a\ for the stochastic process is then given by 5 n n _i 
acting on site i, while the annihilation operator cij is represented by n5 n ^ n+ i (here <5„ jm is the 
Kronecker symbol). In this notation one has \F(n, t)) = J2 n F{n,t)\n). From the definition of the 
time evolution ( |3.1| ) we find that a state vector at time t is given in terms of the vector at time 
t = by \F(n,t)) = exp(-Ht) \F(n,0)). Time-dependent operators A(t) (in euclidean time) are 
defined by A(t) = exp (Ht)Aexp (—Ht). Note that in one dimension this model is equivalent to 



an interface growth model with unrestricted height gradients ||13|| . 

H has a degenerate ground state with energy E = 0, each sector with fixed particle number 
N contains such a vector which are the steady states of the system. Here we shall study vacuum 
expectation values 

i ? (ri,r a ;ti > t2) = <0|a Pl (t 1 )aJ 3 (t2)|0> , (3.3) 

i.e. correlation functions in the vacuum state containing no particles. Here r% t 2 denote sites on the 
ID lattice with lattice constant a. We furthermore assume t\ > £2- In the language of stochastic 
processes this is the conditional probability of finding the particle on site r\ at time t\ provided 
it was on site r2 at time £2- Two-point correlation functions in iV-particle steady states can be 



reduced to such correlators 1114 



In order to compute the correlation function ( |3.3|) we first note (a) that since we are considering 
a steady state correlator, F depends only on rj, T2 and r = t± — ti and (b) that because of 
translational invariance, F depends only r = t\ — r2 and r. Furthermore, operators a\ and a« 
satisfy the equations of motion 

d t a\ = [H,a\] = -\ (a\ +1 + 4-i - 2aJ) ^ 
d t di = [H , a,i] = \ (a i+ x + ai_i - 2aj) . 

This is equivalent to Ca\ = Ca{ = where C is the lattice Schrodinger operator ( |2.6| ) with the 
eigenvalue A4 — 1 of the mass operator M defined by [N , a{] = —Aiai = —ai and N , a\ 
Ma\ = a\. 
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Eqs. ( |3.4| ) are easy to integrate and one finds 

F(r, r) ee F(n, r 2 ; t x , t 2 ) = e^r) (3.5) 



which is in full agreement with the prediction (|2.23|) with amplitude g = 1, mass M = 1, lattice 



constant a = 1 and critical exponent x — 1/2 obtained from Schrodinger invariance and current 
conservation which is implied in eqs. ( |3.4|) . 

3.2 Driven non-exclusive lattice diffusion 

Symmetric lattice diffusion as discussed above appears to be simplest many particle system which 
is invariant under the discretized Schrodinger symmetry introduced in section 2. One may also 
study its driven version where particles hop with different rates pr = qrii and pl = q~ x ni from site 
i to their right and left nearest neighbouring sites i ± 1 respectively. In this case the Hamiltonian 
reads 

H = - {W+iK+i - a d + ^(4+1 - 4) a i+i} (3-6) 

i 

and the correlation function ( |3.3| ) is found by integrating the equations of motion for and a\ to 
be 

F{r, t) = F(r u r 2 ; t u t 2 ) = e"^ e~ T J r (r) (3.7) 

with v = — \nq and 6 = 2 sinh 2 (i;/2) = |(g + q^ 1 ) — 1. The same result can be obtained from yet 
another lattice realization of the Schrodinger algebra. To see this, we go back to the algebra A 
(|2.1| ). Rather than reusing the realization given in eqs. (|2.3| , |2.4|) , we now take 



E = d t + Mb , T = t , P=-smh(^-(d r + Mv)) , X= - 1 - r (3.8) 

a \2 J cosh(|(9 r + Mv)) 

Using the definitions in table 1 this leads to a new realization of the Schrodinger operators and 
one may go again through the calculation of the two-point function. The result is 

F(n, t) = g 6 Ml M^ 1 ^- M ^- Mim t 1 ^e- t I n {t) (3.9) 



with n and t defined as in Q2.24 ). This is reproduced by the result (|3.7|) for asymmetric diffusion 



with go = 1, a = 1, M = 1 and x — 1/2. Note that in the continuum limit a —>■ the correlation 
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function is proportional to exp (— M.v 2 t/2 — M.vp) -t~ 1 / 2 exp {—M.p 2 j (2t)). This is the correlation 
function of the undriven system, but observed in a frame of reference moving with constant velocity 
v. On the lattice, eq. ( |3.9|) is understood as follows. A Galilei transformation will bring us to the 
frame where particles in the mean are at rest. The effect of this Galilei transformation produces 
the result ( |3.9|) , and the condition b = 2a~ 2 sinh 2 (at> /2) takes care that the stationary state has 
indeed vanishing energy. 

3.3 Exclusive symmetric lattice diffusion 

We come back to lattice diffusion as defined in section 3.1 but demand that on each site there can 



be at most one particle [y| . Then the same expression (|3.5| ) for the two-point correlation function 
(|3.3| ) is found. The time evolution operator of this system is given by the spin 1/2 isotropic 
Heisenberg ferromagnet 



H 



7 E (>*X + i + ^W+i + « + i - l) (3-10) 



where a^' y,z are the Pauli matrices. This is an interesting example of a process where the time 
evolution operator is given in terms of operators satisfying a Hecke algebra ||I6| . The equations of 
motion for the creation and annihiliation operators a ± = \{cr x ± ia y ), respectively, are nonlinear 

<V+ = ~ «! + °U - 2af) + (a- +1 a+ +1 + a^a+j) af (3.11) 

and a similar expression for a~ . However, the nonlinear piece is proportional to the number 
operators at sites i ± 1 and therefore vanishes in the expession ( |3.3| ). The linear piece is again the 
lattice Laplacian which accounts for the correlation functions being given by the same expressions 
as in the non-exclusive case. 

We only showed here that the two-point function can be reproduced in the models studied and 
it is not immediately obvious that the same should be true for the higher correlators. But one 
might justify the consideration of Schrodinger invariance in this context by noting that Schrodinger 
invariance can be shown to hold also for at least certain types of non-linear Schrodinger equations, 
e.g. (id t + d 2 + u\ip\ 2 )ip = 0, see |17|, for more information. The equations of motion (|3.11| ) 



are a discretized version of this equation. Furthermore, known results for a partially integrated 
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four-point function of the exclusion process ]T3[, suggest that one might predict such correlators 
from an invariance of the kind discussed here, but adapted to systems with boundaries. 



3.4 Other stochastic models 

The same type of correlation functions can be reproduced from other systems as well. Consider the 
ID kinetic Ising model with Glauber dynamics [19]. The connected spin-spin correlation function 



< s ri (ti)s r2 (t2) > c of spins at the lattice sites n t 2 at times ti,2, when the system is in equilibrium 
at a temperature T, can be calculated exactly. If in the initial state there are no correlations 
between the spins, one has fl9| 



< s ri (t!)s r2 (t 2 ) > c = e- at I r (a^t) (3.12) 

with r = 77 — r 2 , t = \ti—t 2 \, a is a constant reaction rate, 7 = tanh(2J/ (fcgT)) where J is the Ising 
model coupling constant. For 7 = 1 corresponding to the zero temperature static critical point 
of the one-dimensional Ising model this is obviously in agreement with the prediction eq. ( |2.23| ), 
upon identification of the non-universal parameters. But also the off-critical form of the correlation 
function (7 7^ 1) is in agreement with lattice Schrodinger invariance ( |3.9| ) with v — 0, x — 1/2 
and choosing b and Ai appropriately. This corresponds to the invariance group of a Schrodinger 
equation in a constant potential, proportional to b. 



Another example is provided by the stochastic adsorption-desorption process of dimers [[20 
Consider a ID lattice whose sites may be empty or occupied by at most one particle. A pair of 
empty nearest-neighbour sites may become occupied with rate e (dimer adsorption), and a pair of 
occupied nearest-neighbor sites may become empty with the same rate e (dimer desorption). A 
single particle may hop from an occupied site to an empty nearest neighbor site with a rate h. Let 
n r [t) denote the occupation number of site r at time t. It can be shown that the system evolves 
towards a steady state, characterized by < n T >= 1/2. In the steady state, the density- density 



correlation is on an infinite lattice 20 



< n ri (t)n r2 (0) >= ^e- 2{e+h)t I r {2(h - e)t) (3.13) 



with r = ri — r 2 . This correlation again reproduces lattice Schrodinger invariance (|3.9| ), when 
choosing b appropriately. 
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4 Conclusions 



We have introduced new realizations of the Schrodinger algebra as describing dynamical symme- 
tries of discrete lattice systems. The two-point correlation function is found from the requirement 
that it covariantly transforms under these realizations and includes the previously known result of 
the continuum as the special limiting case a — > 0. We find that the two-point function so obtained 
can be reproduced from the two-particle correlation function, calculated in the stationary state, 
of several simple stochastic models. 

The construction of lattice realizations of continuous space-time symmetries using the algebra 
( |2.1| ), ( j2.2|) extends beyond the few cases studied here. Likewise, one may generalize this approach 
to higher dimensions. This will be reported elsewhere. 
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